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Abstract —Multiplication is one of the most important oper- 
ation in computer arithmetic. Many integer operations such as 
squaring, division and computing reciprocal require same order 
of time as multiplication whereas some other operations such 
as computing GCD and residue operation require at most a 
factor of logn time more than multiplication. We propose an 
integer multiplication algorithm using Nikliilam method of Vedic 
mathematics which can be used to multiply two binary numbers 
efficiently. 

Index Terttis —Integer Multiplication, Algoritlim, Computer 
Arithmetic, Vedic Mathematics, Computation 

I. Introduction 

The classical method of adding two integers of n-bits takes 
0(n) bit operations but the classical method of multiplying 
them takes 0(n 2 ) bit operations. Complexity of addition is 
optimal in number of bit operations, whereas optimal multipli- 
cation algorithm for integers is an open problem. The classical 
approach to multiply two n-bit integers requires 0(n 2 ) bit op- 
erations. Karatsuba multiplication [9] uses divide-and-conquer 
technique to multiply two n-bit integers in O(n log3 ) bit oper- 
ations(logarithms are to the base 2 unless otherwise specified) 
by replacing some of the multiplication by less costly addition 
and subtraction. Toom-Cook algorithm further improves the 
above bound [3]. Toom-Cook method is the generalization of 
Karatsuba method which split each number to be multiplied 
in multiple parts [12]. Given two large integers, Toom-Cook 
splits up multiplicand and multiplier into k smaller parts each 
of length l, and performs operations on the parts. As k grows, 
one may combine many of the multiplication sub-operations, 
thus reducing the overall complexity of the algorithm. For k = 
3 Toom-Cook reduces 9 multiplication to 5, with asymptotic 
complexity of O(n( log5 )/( log3 )). Schonhage-Strassen integer 
multiplication algorithm [10] uses Fast Fourier Transform 
(FFT) by selecting the principal roots of unity as evaluation 
point to perform multiplication in 0(nlog?rloglogn) bit op- 
erations. FFT method employs a divide-and-conquer strategy 
by taking advantage of the special properties of the complex 
root of unity to perform multiplication of two polynomials 
in only O(nlogn) arithmetic complexity [4]. In fact many 
multiplication algorithms can viewed as schemes for the 
evaluation of polynomials then multiplication of their values 
and followed by interpolation. Currently, the asymptotically 
fastest algorithm for multiplication of two /t-bit integers is by 
Furer [7] which runs in O(nlogn2 0 ( log O). Where log* n is 
iterated logarithm function [5] defined as: 

log* n = min{« > 0 : log ( '*' ) n < 1}. 


Furer algorithm uses arithmetic over complex number. Same 
asymptotic bound can also be achieved using modular arith- 
metic [ 6 ]. 

The Schonhage-Strassen algorithm and Furer algorithm are 
asymptotically fast but they are suited for extremely large 
numbers. Furer algorithm although asymptotically fastest, only 
achieves an advantage for astronomically large values and 
as such it is currently not used in practice. The crossover 
points between these algorithms are usually very high when 
the algorithms are implemented [ 8 ]. For small inputs even 
Karatsuba algorithm runs slower than the classical multipli- 
cation algorithm because of recursion overhead. In this paper 
we use Nikhilam sutra or method from vedic mathematics [11] 
to perform efficient multiplication for small inputs. Nikhilam 
sutra performs large multiplication by converting it to small 
multiplication along with some addition and shifting opera- 
tions. 

This paper is organized as follows. Section 2 describes 
background and motivation. Section 3 presents our proposed 
work, the multiplication algorithm and its features. In section 
4, we present some applications. Finally, section 5 contains 
conclusion. 

II. Background and Motivation 

The simplest method to multiply two n-digit integers is 
using classical or long multiplication method which requires 
0(n 2 ) multiplication operations. Whereas to add or subtract 
two n-digit integers using traditional method requires at most 
n number of addition or subtraction which is optimal in terms 
of number of addition/subtraction operation performed. To 
improve the 0(n 2 ) bound of multiplication, several algorithms 
have been discovered. The simplest one is Karatsuba algorithm 
which is based on divide-and-conquer paradigm [1]. Karatsuba 
algorithm is based on the fact that two-digit multiplication can 
be done with only three instead of four multiplications required 
by standard method. Suppose we want to multiply two 2-digit 
decimal numbers aia 2 * b\b 2 : 

1. Compute A = ai * bi 

2. Compute B = a-i * 62 

3. Compute C = (ai + 02 ) * (foi + b 2 ) 

4. Compute D = C — A — B, here D is equal to ai*b 2 + 
a 2 * b\ 

5. Result 100 *A +10 *D + B 

For large number of digits we can apply this method 
recursively by splitting the multiplicand and multiplier in two 
parts. The complexity of this method is 0(n^ log3 ^^ log2 ^). 



Since multiplication operation is costly as compared to 
addition and shift, some constant number of addition and shift 
operations can be safely ignored. In this paper we assume that 
multiplicand and multiplier are having equal number of digits. 

For example suppose we want to multiply 95 * 96. The 
standard method of long multiplication requires 4 one-digit 
multiplication along with some addition and shift. Using 
Karatsuba algorithm we can compute it as follows: 

1. Compute A = 9 * 9 

2. Compute B = 5*6 

3. Compute C = (9 + 5) * (9 + 6); C = 14 * 15 

4. Compute D = C - A - B = 210 - 81 - 30 = 99 

5. Result 100 * 81 + 10 * 99 + 30 = 9120 

A schematic view of above multiplication is shown in Fig. 
1. The three multiplication operations are enclosed in ellipse. 
In fact the total number of 1-digit multiplication required in 
above example is 5. Note that to compute 14 * 15 in step 3 
requires, three 1-digit multiplication by applying this method 
one’s more. 



Fig. 1. Multiplication of integers (95 * 96) using Karatsuba method 

Nikhilam Sutra is one of the 16 sutras of Vedic mathematics 
[11]. It can be used to convert large-digits multiplication to 
small-digits multiplication with the help of few extra add, 
subtract and shift operations. In some cases two-digit multipli- 
cation can be performed using only 1 one-digit multiplication 
instead of 3 one-digit multiplication as required by Karatsuba 
algorithm. Suppose we have to perform same multiplication 
95 * 96 using this method. We can use the Nikhilam sutra as 
follows: 

1. Compute A = 100 — 95; Subtract the multiplicand from 
nearest base 

2. Compute B = 100 — 96; Subtract the multiplier from 
the same base 

3. Compute C = B * A = 5*4 = 20 

4. Compute D = 95 - 4 = 96 - 5 = 91 

5. Result 100 * D + C = 9120 

In the Fig. 2, we can see that there is only one multiplication 
operation involved. 

Above multiplication is also shown in Table 1. In this 
multiplication we have used 1 multiplication, 1 addition, 3 
subtraction and 1 shift operation. This particular multiplica- 
tion is more efficient than both standard multiplication and 



Fig. 2. Multiplication of integers (95 * 96) using Nikhilam method 

Karatsuba method. Suppose multiplicand is rn = x — a and 
multiplier is n = x — b where x is nearest base. We have: 

m * n = (x — a) * (x — b) = x(x — a — b) + ab 

The general scheme of multiplication m*n is shown in Table 

2 . 

TABLE I 

Multiplication of 95 * 96 



TABLE II 

Multiplication of m* n 



This scheme can also be utilized, when the multiplicand and 
multiplier are just above the certain power of base(which is 
10 in this case). Suppose we have to multiply 105 * 106. We 
can proceed as follows: 

1. Compute A = 105—100; Subtract the multiplicand from 
nearest base 

2. Compute B = 106 — 100; Subtract the multiplier from 
the same base 

3. Compute C = B * A = 5*6 = 30 

4. Compute D = 105 + 6 = 106 + 5 = 111 

5. Result 100 *D + C = 11130 

Above multiplication is shown in Table 3. This 3-digit 
multiplication is performed in just single 1 digit multiplication. 











TABLE III 

Multiplication of 105 * 106 



Integer 

Base Difference 

Multiplicand 

105 

(105-1003=5 

Multiplier 

106 

(106-1003=6 


(105+6)=111 

(5*6)=30 

Result 

11130 



TABLE V 

Binary Multiplication OF 11 * 11 



Bits 

Base Difference 

Multiplicand 

11 

(11 103=1 

Multiplier 

11 

(11-103=1 


(11+13=100 

(1*1)=1 

Result 

1001 



whereas the standard method will take 9 multiplication and 
Karatsuba can compute in 4 multiplication operations. 

The principle behind this is as follows: Let the multiplicand 
be m = x + a and multiplier be n = x + b where x is the 
nearest base, then 

m * n = (x + a) * (x + b) = x(x + a + b) + ab 

Details are given in Table 4. 

TABLE IV 

Multiplication of m*n 



Integer 

Base Difference 

Multiplicand 

m 

(x + a) — x 

Multiplier 

n 

(x + b) — x 


(x + a + b) 

ab 

Result 

x(x + a + b) + ab 



III. Proposed Work 
A Binary Multiplication 

We can perform binary digit multiplication using Nikbilam 
sutra by converting n-bit multiplication to (n — l)-bit 
multiplication and some additional add/subtract and shift 
operation. We can apply this conversion repeatedly until we 
get trivial multiplicand/multiplier or 1-bit multiplication. We 
can also put some threshold limit m where 1 < m < n up to 
which we would like to do this conversion. 

2-bit multiplication can be performed using single 1-bit 
multiplication. For example if we have to multiply 11 * 11. 
Here multiplicand M = 11, and multiplier N = 11. We can 
proceed as follows: 

1. Compute A = 11 — 10; Subtract the multiplicand from 
nearest base 

2. Compute B = 11 — 10; Subtract the multiplier from the 
same base 

3. Compute C = B * A = 1*1 = 1 

4. Compute D = M + B = N + A= 11 + 1 = 100 

5. Result 10 * D + C = 1001 

The only multiplication required in this computation(Table 5) 
is for C in step 3. 

For 3-bit multiplication consider the example of 101 * 110 


1. Compute A = 101 — 100; Subtract the multiplicand from 
nearest base 

2. Compute B = 110 — 100; Subtract the multiplier from 
the same base 

3. Compute C = 10 * 1 = 10 

4. Compute D = M + B = N + A = 101 + 10 = 111 

5. Result 100 * D + C = 11110 

In this computation also two 1-bit multiplication is performed. 
While in case of standard multiplication 9 multiplication is 
required, and Karatsuba algorithm use 4 multiplication. 

TABLE VI 

3-bit Binary Multiplication of 101 * 110 



Bits 

Base Difference 

Multiplicand 

101 

(101-100)=1 

Multiplier 

110 

(110-1003=10 


(101+10)=111 

(1*10)=10 

Result 

11110 



Let us consider 4-bit multiplication of 1111 * 1111. In this 
case multiplicand(M) and multiplier(N) are equal i.e. M = 
N = 1111. For multiplication we can proceed as follows: 

1. Compute A = 1111 — 1000; Subtract the multiplicand 
from nearest base 

2. Compute B = 1111 — 1000; Subtract the multiplier from 
the same base 

3. Compute C = A — 100 = 11 

4. Compute D = B — 100 = 11 

5. Compute E = C — 10 = 1 

6. Compute F = D — 10 = 1 

7. Compute G = E * F = 1 

8. Compute H = (C + F) * 10 + G = 1001 

9. Compute I = (A + D) * 100 + H = 110001 
10. Result J =(M + B)* 1000 + 1 = 11100001 

In this computation the only multiplication involved is in 
step 7 and apart from that all other steps are having either 
add/subtract or shift operation. Details of this multiplication 
are shown in Table 7. 

From these examples we can easily recognize that the 
computation is simple when multiplicand and multiplier both 
are same. 







TABLE VII 

4-bit Binary Multiplication of 1111 * 1111 



Bits 

Base Difference 

Next Difference 

Next Difference 

Multiplicand 

1111 

(1111-10001=111 

(111-1001=11 

(11-101=1 

Multiplier 

1111 

(1111-10001=111 

(111-1001=11 

(11-101=1 





1*1 = 1 




(11 + 1) * 10+ 1 = 1001 




(111 + 11) * 100 + 1001 = 110001 




(1111 + 111) * 1000 + 110001 




Result 

11100001 





B. NikMlam Algoritkm 

Nikhilam sutra seems to have special advantage when the 
multiplicand and multiplier are same. If both multiplicand and 
multiplier are equal then multiplication operation is known as 
squaring. Squaring is considered as special case of multiplica- 
tion. First we write the NikbilamSquaring algorithm which can 
be used to compute square of a binary integer. NikhilamSquar- 
ing algorithm is in turn called by NikhilamMultiplication to 
perform the multiplication. Multiplication and squaring related 
to each other by the following well known formula: 

{x + y ) 2 - (x-y ) 2 
y 4 

Description of NikhilamSquaring is given in Algorithm 1. It 
takes input .4 as a binary number of n-bits and produces square 
of A as its output which can be up to 2n-bits. First while-loop 
is used for the computation of the forward direction subtraction 
operations. Two counters i and j are used to keep track of 
processed input and proper base subtraction respectively. If- 
loop is used to check whether the corresponding bit is 0 or 
1. If base power to be subtracted is more than the number 
itself then else part of if-loop is executed. This happens to be 
only when the first bit of the number is 0. Least significant 
bits multiplication is assigned in B \. Second while-loop is 
used for the computation of reverse direction for shifting and 
addition operations. Again If-loop is used to check whether 
two consecutive Aj values are same, if it is so next value of 
Bi is unchanged, otherwise its value is updated. 

Execution of NikhilamSquaring algorithm for 6-bit input 
A = 101010 is shown in Table 8. We have Ai = 101010, 
A 2 = (101010 - 100000) = 01010, A 3 = 1010 (since 
01010 > 10000), A 4 = (1010-1000) = 010, A 5 = 10, A 6 = 
o, B ± = 0 * 0 = 0, B 2 = (10 + 0) * 10 + 0 = 100, B 3 = 100 
(since A 4 = A 5 ), B 4 = (1010 + 010)* 1000 +100 = 1100100, 
B 5 = 1100100, B 6 = (101010 + 1010)100000 + 1100100 = 
11011100100 . 

Correctness of the NikhilamSquaring can be easily 
established using induction on the bit length of the input 
A. Note that in each column of Table 8, the partial result 
is, in fact multiplication of corresponding multiplicand and 
multiplier in that column. 


Algorithm 1 : NikhilamSquaring (A) 

INPUT: A = £o -1 o.,x l 
OUTPUT: B = A*A = £o" _1 b k x k 
A x <- A 

i <— 2, j ^— ti 1 

while (i <n and j > 1) do 

if (A t > 2>) then 
Ai <— Ai -1 — 2 J 

else 

Ai <— i 

end if 

i <- i + l,j <- j - 1 

end while 

B\ = A n * A n 
i <— 2, j <— n — 1 

while (i <n and j > 1) do 
if Aj / Aj o -| then 

Bi <— Bi-i + ( Aj + Aj+i)2 l 1 

else 

B, <— Bi_ i 

end if 

i <- i + l,j <- j - 1 

end while 
return B B n 


Theorem 1. NikMlamSquaring algoritMn computes square 
of the input A. 

Proof: We prove this using induction on the bit length n of 
A. For n = 1, number A has only one bit and ,4 2 1 when 

A is 1 and A 2 = 0 when A is 0. Therefore it works for n = 1. 
Assume it is true for k — 1. Now we show it for k. Assume 
A k _i = x, therefore A 2 _^ = x 2 . 

Case 1: If kth bit is 1 then A k = 2x + 1, and the processing 
of fcth step is 

((2x + 1) + x) * 2 k + x 2 

also we have 


(2x + 1) - 2 k = x => 2 k = x + 1 




TABLE VIII 

6-bit Binary Multiplication of 101010 * 101010 


Binary Digits Base Difference Next Difference Next Difference Next Difference Next Difference 


Multiplicand 

A\ = 101010 

A 2 = 01010 

A 3 = 1010 

A 4 = 010 

o 

II 

lO 

Ae =0 

Multiplier 

A\ = 101010 

A 2 = 01010 

A 3 = 1010 

A 4 = 010 

o 

II 

to 

Ag = 0 

B\ = Aq * Aq 






B 4 = 0 * 0 = 0 

b 2 





B 2 = 100 


b 3 




B 3 = 100 



Bi 



B 4 = 1100100 




b 5 


B 5 = 1100100 





Result 

Bq = 11011100100 







and therefore 

((2x + 1) + x) * 2 k + x 2 = ((2x + 1) + x)(x + 1) + x 2 

=> (3x + l)(a; + 1) + x 2 = 4x 2 + 4a; + 1 = (2x + l) 2 

and hence Af = (2x + l) 2 , and the statement of theorem 
follows. 

Case 2: If fcth bit is 0, then the statement is trivial and 
Ak = = x 2 , and theorem is proved. 

NikhilamMultiplication is described in Algorithm 2. It 
takes two binary numbers A and B as input and compute 
their multiplication(C') as output. It performs two calls 
to NikhilamSquaring algorithm and using that, outputs 
the desired result. NikhilamMultiplication uses only 1 
multiplication, 1 division and 0(c* n) add/subtract and shift 
operations for some constant c. To show the correctness of 
NikhilamMultiplication is trivial given the correctness of 
NikhilamSquaring. 

Algorithm 2 : NikhilaniMultiplication(,4, B) 

INPUT: A = X)o _1 a i xi > B = Eo” 1 b i xi 
OUTPUT: C = AB = Eo”” 1 c kX k 

Di c- NikhilamSquaring (A + B) 
l)‘2 <r- NikhilamSquaring (A — B) 
return C <= (D\ — D 2 )/A 

Theorem 2. NikhilamMultiplication algorithm computes 
the product of A and B. 

Proof: The statement of the above theorem follows from 
Theorem 1 and the fact that : 

A, B = ( A + B ) a -i A - B ) a 
4 


IV. Applications 

Since, asymptotically fast multiplication algorithms like 
Schonhage-Strassen and Furer algorithms are only useful for 
extremely large numbers, for small to medium size numbers 
we can apply Nikhilam multiplication. We can use Nikhilam 
multiplication even in conjunction with some other fast algo- 
rithm like Karatsuba. If n 0 is the threshold between classical 
multiplication and Karatsuba algoritfun then up to threshold 
limit n 0 we can use Nikbilam multiply and beyond that limit 
we can use Karatsuba multiply. We can write the Karatsuba 
algorithm as given in [2], with the only modification that if 
n < Uq NikhilamMultiplication is called. The corresponding 
pseudo code is given in Algorithm 3. Karatsuba multiplication 
has relatively small threshold with the classical multiplication. 
The optimal threshold for Karatsuba algorithm can vary from 
about ten to hundred words. NikhilamMultiplication can also 
be used as a stand alone multiplication algorithm. 

Further optimization to NikfulamMultiplication is also pos- 
sible. Least significant 0’s can be truncated from input num- 
bers, and in the end corresponding modification can be done 
in single operation. We can process the multiplicand and 
multiplier for consecutive 0’s to skip some of the steps of 
the algorithm. The proposed algorithm is particularly efficient 
because multiplication operation is least involved in it. 

Algorithm 3 : KaratsubaMultiplication(A, B) 

INPUT: A = Eo” 1 ctix\ B = E 0 _1 b o xi 
OUTPUT: C = AB = Eo"” 1 c kX k 

if (n < n 0 ) then 

return NikhilamMultiplication (A,B) 
end if 

K <= (n/ 2 ) 

(Aq, Bo) <— (A, B) mod x k , (A\,B\) <- (A, B) div x k 
sa sign(A 0 - Ai), s B <= sign (B 0 - B\) 

C 0 <- KaratsubaMultiplication (A 0 ,B 0 ) 

C\ <- KaratsubaMultiplication (A\,B\) 

C 2 <- KaratsubaMultiplication (|/1 (J A\\.\B tl B/) 

return C <— C 0 + (C 0 + C\ — SAS B C 2 )x k + C\x 2k 




V. CONCLUSIONS AND FUTURE WORK 

In this paper we have explored the possibility of applying 
the Nikhilam sutra of Vedic mathematics to binary number 
multiplication. We can take advantage of the fact that this sutra 
can convert large-digit multiplication to corresponding small 
digit multiplication. Nikhilam method is particularly efficient 
when both multiplicand and multiplier are near to some base 
(radix) power. To take this advantage, we have first performed 
square operation in NikhilamSquaring and then we have used 
this to finally compute multiplication. Future work can be 
to extend this method to large-digit multiplication and exploit 
it’s properties to perform fast integer multiplications. 
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